On asymptotic stability of the Skyrmion 
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We study the asymptotic behavior of spherically symmetric solutions in the Skyrme model. We 
show that the relaxation to the degree-one soliton (called the Skyrmion) has a universal form of 
a superposition of two effects: exponentially damped oscillations (the quasinormal ringing) and a 
power law decay (the tail) . The quasinormal ringing, which dominates the dynamics for intermediate 
times, is a linear resonance effect. In contrast, the polynomial tail, which becomes uncovered at late 
times, is shown to be a nonlinear phenomenon. 



Introduction. Stable stationary solutions are natural 
candidates for the endstates of evolution of many phys- 
ical systems. The relaxation to these states is well un- 
derstood for dissipative systems, described by diffusion 
equations, however for conservative Hamiltonian systems 
on unbounded domains the problem is much more diffi- 
cult because there is no local dissipation of energy and 
convergence to equilibrium is due to dispersion, that is, 
radiation of excess energy to infinity Understanding 
dissipation by dispersion is important physically because 
the radiation emitted during the approach to equilibrium 
encodes information about the attractor - this kind of 
inverse problem has various applications, for instance in 
identifying black holes via gravitation radiation emitted 
during the last stages of gravitational collapse. 

In this paper we address the problem of relaxation to 
equilibrium in a very simple setting of the spherically 
symmetric Skyrme model 2] . This model, apart from its 
physical relevance in particle physics, is attractive theo- 
retically because the equilibrium state is completely rigid: 
it has no moduli and no internal degrees of freedom, 
which makes the mathematical analysis feasible. Yet, de- 
spite the simplicity of the model, the relaxation process 
exhibits a surprising feature: after a transient oscillatory 
exponential decay (so called quasinormal ringing) which 
is a linear resonance effect, there proceeds a power law 
tail which has a nonlinear origin. Pointing out the failure 
of the linear perturbation theory in capturing the asymp- 
totic dynamics is the main message of this paper. Below, 
after introducing the model, we first describe the quasi- 
normal ringing using the linear perturbation theory, then 
we present the numerical evidence for the asymptotic be- 
havior of solutions, and finally we explain the tail using 
the nonlinear perturbation theory. 

Background. Let M be a spacetime with a metric r/^ 
and N be a complete Riemannian manifold with a met- 
ric gAB- Consider a map U : M — » N and denote by 
Sfj,u = gAB dfJJ A d v U B the pulled back metric. The (gen- 
eralized) Skyrme model is defined by the lagrangian 



where a is the coupling constant having the dimension 
of length. In this paper we consider the original Skyrme 
model 0] where M is the 3 + 1 dimensional Minkowski 
spacetime with the metric r\ = — dt 2 + dr 2 + r 2 dio 2 
and N is the 3-sphere with the round metric ds 2 — 
dF 2 + sin 2 FdQ 2 , where dio 2 and dQ 2 are the standard 
metrics on the unit 2-sphere. We restrict our attention 
to corotational maps for which F = F(t,r) and fi = u>. 
For such maps the Euler-Lagrange equations correspond- 
ing to |T]) reduce to the single nonlinear wave equation 
(using the abbreviation w = r 2 + 2a 2 sin 2 F) 



(wF)-(wF')' +$m(2F)+a 2 sin(2F) 



sin 2 i^ 



(2) 

where primes and dots denote derivatives with respect to 
r and t, respectively. We are interested in the long time 
behavior of solutions of this equation for smooth finite 
energy initial data. Regularity at the origin is ensured by 
the boundary condition F(t, 0) = 0. The total conserved 
energy associated with solutions of equation ([2]) can be 
written as the sum E = E a + E$, where 
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The quadratic part of the energy E a , corresponding to 
the pure sigma model (a — 0), has the supercritical 
scaling E a [F{x/\)\ = XE a [F(x)], hence for a = it 
may be energetically favorable for solutions to shrink 
and consequently singularities are expected to develop 
for some initial data [4| . The quartic part of the energy 
Es, introduced by Skyrme, has the subcritical scaling 
Es[F(x/X)] = \~ 1 E s [F(x)], thus for nonzero a shrink- 
ing of solutions to zero size is prevented by energy con- 
servation. Hereafter, we assume that the coupling con- 
stant a is nonzero and use the unit of length such that 
a = 1. Note that for the total energy to be finite, so- 
lutions must satisfy the boundary condition at spatial 
infinity F(t, oo) = rmr (m = 0, 1, . . . ), where an integer 
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Tji is the topological degree of the map. Since the time 
evolution is continuous (as long as no singularity forms), 
this condition breaks the initial value problem into in- 
finitely many disjoint topological sectors labeled by the 
degree 77i. It is well known that for each m there is a 
unique regular static solution of equation ^ [B| . Our nu- 
merical studies indicate that these static solutions play 
the role of global attractors in the evolution of regular 
corotational initial data of a given degree, that is every 
solution starting from smooth finite energy initial data 
of degree m remains globally regular for all times and 
asymptotically settles down to the static solution of de- 
gree m. Below, for concreteness, we focus our attention 
on the degree-one sector m = 1. In this case we shall refer 
to the static solution as the Skyrmion and denote it by 
S(r). The Skyrmion is the most interesting corotational 
soliton because, in contrast to solutions with m > 1, it 
is stable with respect to general (nonradial) perturba- 
tions. The profile function S(r) (see Fig. 1) is not known 
explicitly. The existence of the Skyrmion was proved rig- 
orously both by ODE techniques || and by variational 
methods @. 

Linear stability and quasinormal modes. In order to 
interpret the numerical results shown below we first need 
to discuss the spectrum the linear perturbations around 
the Skyrmion. To this end we seek solutions in the form 
F(t,r) = S(r) + 5F(t,r). Plugging this into equation J2), 
linearizing and using the auxiliary field v(t, r) defined by 

v(t, r) 



5F(t,r)= _. 

\l v + 2 sin 2 S 

we get the linear wave equation for the perturbation 
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where the effective potential is 
, 1 + 3a 2 + 3a 4 
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Near the origin S(r) ~ br (with b ks 2.0075) and near 
spatial infinity S(r) ~ ir — c/r 2 (with c w 2.1596), hence 
the potential V(r) is finite at r = and falls off as r~ 6 
for large r. Substituting v(t, r) — e~ i/)(r) into equation 
([6|) we obtain the / = 1 radial Schrodinger equation 



ip" + ( -j + V(r) J V = k 2 i]j . 



(8) 



It is known that this equation has no bound states (and 
the spectrum is purely continuous k 2 > 0) which implies 
that the Skyrmion is linearly stable [7|. 

Although linear stability is an important property of 
a soliton, it provides little information about the asymp- 
totic behavior of solutions near the soliton. The key con- 
cept in the studies of the asymptotic stability of the soli- 
ton is the notion of the quasinormal mode. The quasi- 
normal mode (known also as the resonance) is a regular 



solution of equation ((SJ which satisfies the outgoing wave 
condition for r — > oo 



V>(r) 
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k = n-ir, r > o. 
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The quasinormal mode with the least damping factor T 
is expected to dominate an intermediate stage of the re- 
laxation to the soliton. In order to find this mode we 
use a shooting method which goes as follows. First, we 
express ip in the amplitude-phase form as tp = Aexp(i0) 
and rewrite equation (J8]) as the following system 



A" + A(j>' 2 + 



V + T 2 - ft 2 ) A = 



(10a) 



Ac))" + 2AV - 2QTA = . (10b) 
To ensure regularity at the center we require that 



A(r) 



and 



(r) - — r 



for r -> 0. (11) 



We need to find and T such that the condition ([9]) is sat- 
isfied. A naive shooting method does not work because 
an unwanted ingoing wave contamination of the condi- 
tion ([9} decreases exponentially with r and cannot be 
tracked numerically. To overcome this difficulty (which 
is intrinsic to the problem and any numerical method has 
to cope with it) we first solve equations (|10a|) and (jlObj) 
with the initial condition (|11|) up to some relatively small 
intermediate ro- Next, for r > r$ we define the logarith- 
mic derivative g = ip' /if) = A' /A + i<f>' and replace the 
Schrodinger equation ((SJ by the Ricatti equation 
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We solve this equation backwards in r from some large 
R to ro starting with the initial value 



g{R) 



kh[(kR) 
h^kR) ' 



(13) 



where the Ricatti- Hankel function h\{kr) = {— i+ -g : )e lkr 
is the exact outgoing wave solution of the free (V — 0) Ri- 
catti equation ([12]). The value of R should be chosen suf- 
ficiently large so that the influence of the potential V(R) 
is negligible, however in practice R should not be too 
large in order to avoid numerical instabilities. Matching 
the logarithmic derivatives at the midpoint ro we found 
the quasinormal mode at k = 0.61 — 0.26z. 

Numerics. We solved equation @ numerically for dif- 
ferent degree-one initial data and found that all solutions 
remain globally regular and asymptotically settle down to 
the Skyrmion (see Fig. 1). Proving this asymptotic com- 
pleteness property is a challenging open problem which 
we do not pursue here but take it as the starting point 
for further discussion. 
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FIG. 1: A solution of degree one (solid line) is shown to con- 
verge to the Skyrmion (dashed line). 



Our aim is to understand the asymptotic dynamics of 
convergence to the Skyrmion. Here we mean convergence 
in the pointwise sense, that is we consider the behavior 
of a solution F (t, r) for large t at a fixed distance r. The 
basic mechanism of decay is, of course, dispersion - it 
is clearly seen in Fig. 1 how the excess energy is being 
radiated away to infinity as the solution approaches the 
Skyrmion. During the relaxation process one can distin- 
guish two universal stages of evolution: the oscillatory 
exponential decay which we shall refer to as the quasi- 
normal ringing and the polynomial decay which we shall 
refer to as the tail. The quasinormal ringing is a well- 
known effect in the linear scattering theory [a Q . As the 
name indicates, it is due to the presence of a quasinormal 
mode. The numerical confirmation of this fact is given 
in Fig. 2 where we show that for intermediate times the 
deviation of the solution from the Skyrmion is perfectly 
approximated by the fundamental quasinormal mode. 

For later times the quasinormal mode becomes negligi- 
ble and the decay takes the form of a polynomial tail. The 
numerical computation of the tail (see Fig. 3) gives the 
power-law decay F(t,r) — S(r) ~ t~ 5 . Surprisingly, the 
rate of decay of the tail is different than that predicted 
by the linear scattering theory. To see this, recall that 
according to this theory [l(J, for compactly supported 
initial data a solution of the linear radial wave equation 
with a regular potential V(r), 

v-v"+ ^ + Vv = , (14) 

decays at a fixed r as v(t,r) ~ £~ 7 , where 7 = 21 + (3 
and (3 > 3 is the rate of fall-off of the potential at spa- 
tial infinity, i.e. V(r) ~ r~@ for r — ► 00 (for compactly 
supported or exponentially localized potentials there are 
no tails). In the case at hand (equation ((6])) we have 
/ = 1 and (3 — 6, hence according to the linear theory 
the tail should have the power index 7 = 8 instead of 
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FIG. 2: We plot ln\P(t,r = 10) ] , where P(t,r) = wF(t,r). 
Fitting the exponentially damped oscillation P(t,ro) = 
Ae~ Ft s'm(Slt + S) to the numerical data on the time interval 
(20, 60) we get SI = 0.610 and F = 0.260 in perfect agree- 
ment with the perturbative calculation of the fundamental 
quasinormal mode. 



5 




-35 1 ' 1 ' — ' 1 1 — 3 

10 100 1000 

t 



FIG. 3: The same plot as in Fig. 2 but using the log-log 
scale. Fitting the function a — bint + c/t (the power law 
decay plus the first correction) to the numerical data on the 
interval (200, 1000) we get b = 6.05. 



the observed 7 = 5. Thus, the tail must be a nonlinear 
effect and to understand it one needs to go beyond the 
linear perturbation theory. We found that the third order 
perturbation theory provides a very good approximation 
of the tail. Unfortunately, the perturbative expansion 
around the Skyrmion is very messy and the technical de- 
tails of the calculation might obscure the key mechanism 
which is responsible for the nonlinear tail. To avoid that, 
we will take advantage of the fact that the same mech- 
anism is operating in a simpler setting of the relaxation 
to the vacuum in the topologically trivial sector m = 0. 
In this case there is no tail at all at the linear level. 
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Nonlinear tail. For topologically trivial initial data all 
solutions converge asymptotically to the vacuum Fq = 0. 
To determine the rate of convergence, we substitute the 
expansion F = eF 1 +e 2 F 2 +e 3 F 3 +0(e 4 ) into equation ©, 
where eF\ satisfies initial data while all F n with n > 1 
have zero data. In the first order we get the free I = 1 
radial wave equation 



CF X = , 



C := 



9*_ 
dt 2 



i_d_ 

r 2 dr 



,d_ 

dr 



(15) 



whose general regular solution has the form 

= a'(t-r)+a'(t + r) + a(t - r) - a(t + r) 

(16) 

where the function a(r) is determined by initial data. We 
assume that the initial data have compact support, hence 
F-y has no tail in agreement with Huygens' principle. 

In the second order CF 2 = 0, hence F 2 vanishes, but 
in the third order we get the inhomogeneous equation 



CF 3 



A F 3 

3r 2 1 



(17) 



(18) 



We solve equation ([17]) using the retarded Green's func- 
tion of the operator C 



>\2 



G(t-t',r,r') = [\r-r'\ < t-t' < r+r'] —— } > . 

4r 2 

(19) 

It follows from (fTo| that for large r the term h is of 
lower order in comparison with the first term on the right 
hand side of equation (|17p . thus we can drop it without 
affecting the leading order asymptotics of the tail. Then, 
using double null coordinates u — t' — r' , v = t' + r' , we 
obtain 



t+r 



F 3 (t,r) = 



3r 2 



dv 



(v - t)(t - u) + r 2 



\t-r\ 



(v — u) 2 



Ff(u,v)du. 

(20) 

We are interested in the asymptotic behavior of F^it^r) 
for a fixed r and t — > oo (time-like infinity). Since the 
initial data have compact support, in this limit we can 
change the order of integration in (|2"01) and perform the 
integration over v explicitly. In the leading order we get 



F 3 (t,r)~crr 



64 



a'{ufdu, (21) 



where the constant c is the only trace of initial data. 
We have verified numerically that this formula provides 
a very good approximation of the tail for solutions having 
sufficiently small initial data (see Fig. 4). 



Conclusions. It should be clear from the above discus- 
sion that the nonlinear tail is not an exceptional feature 
of the Skyrme model but it is a general phenomenon in 
scattering theory for nonlinear wave equations which will 
be present whenever the backscattering due to an effec- 
tive potential around the attractor is weaker that the 
backscattering due to a nonlinearity. This kind of phe- 
nomenon does not seem to have been explored in the 
literature and we hope that our letter will initiate in- 
vestigations of an interplay between linear and nonlinear 
effects in relaxation processes. 




FIG. 4: We plot the degree zero solution (solid line) for 
initial data F(0, r) — r z e~ r , F(0, r) — and superimpose 
(dashed line) the analytic prediction for the tail (|21|l with 
c = 35\/37r/1458 0.0737 (note that there is no adjustable 
parameter) . The relative error (representing the contribution 
from higher order iterations) is ~ 5%. 
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